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Abstract: a reversible version of the Navier Stokes equation is studied. A con- 
jecture emerges stating the equivalence between the reversible equation and the 
usual Navier Stokes equation. The latter appears as a statement of ensembles 
equivalence in the limit of infinite Reynolds number, which plays the role of the 
thermodynamic limit. 
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In reference [GCl] a scheme for applying the chaoticity principle, also intro- 
duced there, to fluid motions was discussed, and in [G5] some other simple 
consequences were derived. 

We begin by remarking that the usual derivation of the NS equations assumes 
on an empirical basis that there is a viscous force between fluid layers flowing 
at different speeds, and transforming kinetic energy into heat: it is desirable 
that such an assumption be justified microscopically, [Sp]. On the other hand 
it is difficult to see how a microscopically reversible equation of motion, non 
trivially driven by a non conservative force, could possibly lead to macroscop- 
ically irreversible equations of the NS type, with energy converted into heat in 
a stationary fashion. Unless a dissipation mechanism is a priori assumed at a 
microscopic level. 

One could be tempted to say that irreversible macroscopic motion arises be- 
cause of a mechanism like the one discussed by Grad and Lanford in their 
derivation of the Boltzmann equation. Note however that, although the Boltz- 
mann equation is certainly irreversible, it seems that it cannot be used to derive 
the NS equations in a situation in which there are non conservative external 
forces acting on the system, at least not without extra assumptions besides the 
ones one is willing to assume in deriving the Boltzmann equation itself in a 
system driven to ec^uilibrium.^ 

The physical problem that is stressed here, among many, is that the work 
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Note also that such a derivation will necessarily involve proving an existence theorem for 
the Boltzmann equation of a confined system for a very long (in fact I believe infinite) time 
and this is out of reach of the presently known techniques, [La]. 
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performed by the forces on the system wiU not be dissipated and the system 
will heat up indefinitely, without reaching a stationary state, hence without 
reaching a state described by a stationary probability distribution for the NS 
equation. 

Starting (to fix the ideas) from the reversible microscopic dynamics, the kind 
of "extra assumptions" consist essentially in supposing suitably large ratios 
between the microscopic time and length scales and the macroscopic ones: in 
order to make effectively infinite the time scale over which heat is created. If 
the scale ratios are very large and their relative values also suitably adjusted 
then the motions may be "well approximated" by a dissipative equation and 
reach a stationary state (very slowly drifting away from itself). Strictly speaking 
the equations will still be reversible (involving perhaps also other fields like the 
temperature field) and will be governed by approximately defined transport 
coefficients equal in the average (over time and space) to the observed ones. 

The estimate of the size of the various scale ratios for this picture to be of 
practical interest, and the corresponding evaluation of the errors, is a difficult 
problem that, as far as I know, has not only no clear solution but the methods 
for its solution do not seem to exist yet, [Sp]. 

The above remarks also show that the macroscopic dissipative equations may 
perhaps be equivalent to reversible equations: the aim of this paper is to argue 
about the possibility of notions like "non equilibrium ensembles" and of inter- 
preting the equivalence statement as a statement about ensembles equivalence. 

For the above reasons it will be useful to imagine an idealized dissipation 
mechanism. Here we shall consider a dissipation due to the interaction of 
the system with a "thermostat" which will be modeled by a force acting on 
the system: the function of such forces will be that of imposing a constant 
energy dissipation rate. The force will be determined by Gauss' least constraint 
principle, [LA], [W].^ 

It is easy to derive the equations of a fluid in such a situation. We suppose 
that: 

(1) the fluid is in a periodic box Q with side L, 

(2) that it is incompressible with density p, 

(3) that the energy dissipation is: 



where w = d A u is the vorticity field and is a parameter equal to the 
phenomenological {i.e. experimentally determined) dynamical viscosity: so that 
v = T]p~^ is the viscosity. 

Imposing that the fluid is not viscous but subject to the constraint that (1.1) 
is a constant of motion, implies that the Euler equations for the velocity and 
pressure fields u,p are modified into: 

This type of thcrmostatting force is essentially the well known Nose-Hoover thermostat 
and it has been studied in particular in the fundamental paper [ECM2] which led to the 
formulation of the chaotic hypothesis as a reinterpretation of the earlier Ruelle's principle, 




(1.1) 



[R2], [G2], [GC2]. 
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where A is the Laplace operator and a is not the viscosity but rather it is 
the multiplier necessary to impose the constraint that (1.1) is an integral of 
motion. A simple computation shows that: 

J {d A g ■ w + uj ■ {w ■ du)) dx 

This has odd symmetry in u, so that (1.2) is reversible: if Vt is the flow 
describing the equation solution (so that t ^ Vtu = u{t) is the solution with 
initial data u ) then the transformation i : « — > — u anticommutes with the 
time evolution Vt: 



%Vt = V-t% (1.4) 

The reversible equation (1.2), (1.3) will be called the Gaussian Navicr Stokes 
equation, or GNS. 

A similar equation, with a constraint on the energy contained in each "mo- 
mentum shell" to be constant was considered in [ZJ], which is the first paper 
in which the idea of a reversible Navier Stokes equation is advanced and stud- 
ied. The energy content of each "momentum shell" was fixed to be the value 
predicted by Kolmogorov theory, [LL] . 

Existence of global solutions to the equations (1.2), (1.3) is non trivial even if 
the initial datum Uq is C°°. In view of the conjecture that follows this might 
be a reassuring tribute to the principle of conservation of difficulties. 

In fact let the velocity field u ='^^^0 1 ~ be represented in Fourier 

series with 7 j, = 7^_^. and fc • 7 =0 (incompressibility condition); here k 

has components that are integer multiples of the "lowest momentum" fco = 
Then consider the equation: 



.2 



Efc^o k fffc • 7, 

OLi — 



a = ai + ae, ae = — — 21 ~ 12 ~ ^^'^^ 



where the k 's take only the values < | fc | < K for some momentum cut-off 
K > and 11 ^ is the orthogonal projection on the plane orthogonal to k . This 
is an equation that defines a "truncation on the momentum sphere K of the 
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GNS equations". The actual GNS equations (1.2), (1.3) have the same form 
with K = +00. 

For simplicity we suppose that the mode = is absent, i.e. ^ ^ = Q: this 

can be done if, as we suppose, the external force g does not have a zero mode 
component {i.e. if it has zero average). 

Calling u the solution to the equation (1.5) with initial datum u e C°°, 
and with components with mode k greater than K set equal to zero, it should 
follow that u x^oo' ^— uniformly in x and in t for t in any bounded 
interval. This is however not so easy and in fact I do not know whether (1.5) 
admits a global solution that can be constructed in this way. ^ 

For the purpose of comparison with the NS equation we shall therefore use the 
equations (1.5) and compare them with the corresponding truncation of the NS 
equation (with constant viscosity u) subject to the non conservative field g : 

with the same K. We shall take K very large, fixing it once and for all and 
dropping it from the places where it would belong as a label (to simplify the 
notations). 

In order that the resulting cut-off equations be physically acceptable, and 
supposing that g ^ ^ only for | fc | ~ /co, one shall have to fix K much larger 

than the Kolmogorov scale K = (e;ji/~^)^/"^, where £/j is the average dissipation 
of the solutions to (1.6) with K = +oo, determined on the basis of heuristic 

dimensional considerations by Sk ^ y'] g \^L: see [LL]. 

We shall also call the truncated equations (1.5), (1.6) respectively GNS and 
NS equations omitting the "truncated" : below we always refer to the truncated 
equations, unless otherwise stated. 

The solutions of the equations (1.6) with initial datum u will be denoted 
V^^ u . And we shall admit that there is a unique stationary distribution 
describing the statistics of all initial data u that are randomly chosen with 
a "Liouvillc distribution", i.e. with a distribution jiQid^) proportional to the 

volume measure III fe |<a' '^7 ■ 

This means that given "any observable" F on the phase space T (of the 
velocity fields with momentum cut-off K) it is: 

lim 1 f F{yri)dt= [ F{^')nns{d^') (1.7) 



Examining tiic tiicory of Lcray, [L], for the NS equations one realizes that the "only" diffi- 
culty there was the absence of a uniform a priori estimate on the total vorticity J fii^ dx : in 
the equation (1.2), (1.3) the uniform boundedness of tiie total vorticity is simply automatic 
(as it is a constant of motion). But the theory of Lcray also made essential use of the con- 
stancy of the viscosity coefficient, or more precisely of its positivity and boundedness away 
from 0, see [Gl]. The latter property is false for the coefficient a which, by reversibility, 
takes values of either sign (although, see below, we expect that it is "more often" positive 
than negative on "typical" motions). 
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for all choices of 7 except a set of zero Liouvillc measure. The distribution 
Uns will be called the SRB distribution for the Eq. (1.6) (see the "zero-th law" 
in [UF], [GC2]). 

We make the corresponding assumptions for the solutions of the cut-off GNS 
equations (1.5) (with the same cut-off) and the relative flow Vf"''; we call /jgns 
the associated SRB distribution, [R2], [ER]. 

In particular we shall study the observable o"(7) that we call the entropy 
production rate and that is given by the divergence of the r.h.s. of Eq. (1.5). 

If Dk is the number of modes k with < \k \ < K then the number of 
(independent) components of {7^} is ^Dk and, see (1.5), setting Dk = 

E|fc|<if2i^, one finds: 

a = 2Dk {ai + ag) + 0.1-0.^ = 2DKa + ai-ae (1.8) 
where ai,ae are suitably defined: e.g. 

Cte = — — 2 = = T ^— {1-9) 

E,fc%,P (E,A%J^)^ 

so that a ~ 2DkCx. 

We call cr„s,£„s the time averages of a and of pi^L^J2k 1^1^17 (total 
vorticity, see (1.1)) with respect to the distribution fins- As mentioned we 
expect, on the basis of dimensional analysis, that e„s ~ £fc = yj\ 9 see 
[LL]). 

Correpondingly we consider the solutions of the GNS equations with total 
vorticity Egns = £ = Sns and call (jgns the average of a with respect to jigns- 

The H theorem of RucUe, [R3], tells us that (Jgns > 0, and Ugns > if the 
distribution figns is not proportional to the volume measure. Here we shall 
suppose that this is the case if ,g 7^ : at least if the force field is such that it 
generates a chaotic motion {i.e. a motion with at least one positive Lyapunov 
exponent) for all but a set of zero volume initial data. In this situation: 

Conjecture ("equivalence of dynamical ensembles"): Given e = Sgns suppose 
that the parameter v in the NS equation is adjusted so that Sns — £gns then 
the averages (Tgns , <Jns of a with respect to figns and with respect to /ins and, 
in fact, the averages of all (reasonable) observables coincide, up to corrections 
0{R-^), R = e^/'^L^/^u-\ 

We recognize in this conjecture a statement very analogous to the familiar 
statements on the equivalence of thermodynamic ensembles, with the thermo- 
dynamic limit replaced by the limit i? ^ cxd of infinite Reynolds number. The 
paper [ZJ] contains various statements that can be seen as steps towards the 
formalization given above of the general conjecture. 

It can be weakened very much for the purposes of possible applications: we 
state it in the above very strong form to stress an extreme thought. Applica- 
tions will be envisaged elswhere. 
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On heuristic grounds, it would be justified if one did accept that the entropy 
creation rate reaches its average on a time scale that is fast compared to the 
hydrodynamical scales. The coefficient a. ~ {2DK)~^cr, see (1.8), would be 
confused with its average {a)gns and identified with the viscosity constant v. 

In this way the GNS and the NS equations would be equivalently good: both 
being the macroscopic trace of two equivalent microscopic dissipation mecha- 
nisms: one explicitly specified by the gaussian constraint of constant dissipa- 
tion and the other unspecified but phenomenologically modeled by a constant 
viscosity. 

The interest of the conjecture is that it allows us to deduce properties of the 
"usual" NS equation from properties of the GNS equation: the latter, being 
reversible, can be studied via the chaotic hypothesis of [GCl], [GC2] which 
leads to non trivial predictions like the fluctuation theorem of [GCl]. 
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